EXTINCTION OF FLEMING-VIOT-TYPE PARTICLE SYSTEMS WITH 

STRONG DRIFT 



MARIUSZ BIENIEK, KRZYSZTOF BURDZY AND SOUMIK PAL 

Abstract. We consider a Fleming- Viot- type particle system consisting of indepen- 
dently moving particles that are killed on the boundary of a domain. At the time of 
death of a particle, another particle branches. If there are only two particles and the 
underlying motion is a Bessel process on (0, oo), both particles converge to at a finite 
■ time if and only if the dimension of the Bessel process is less than 0. If the underlying 

diffusion is Brownian motion with a drift stronger than (but arbitrarily close to, in a 



O 

m 



suitable sense) the drift of a Bessel process, all particles converge to at a finite time, 
for any number of particles. 

1. Introduction 

Our paper is motivated by an open problem concerning extinction in a finite time of a 
brandling particle system. We prove two results that are related to the original problem 
and might shed some light on the still unanswered question. 

> , 

OO , The following Fleming- Viot-type particle system was studied in |5J. Consider an open 

O I bounded set D C M°' and an integer N > 2. Let = (X/, . . . ^X^) be a process with 
values in defined as follows. Let Xq = G . Then the processes 

X/, . . . , X^^ evolve as independent Brownian motions until the time ri when one of them, 
say, X^ hits the boundary oi D. At this time one of the remaining particles is chosen 
^ \ uniformly, say, X^ , and the process X^ jumps at time ri to X^^ . The processes X}, . . . , X^ 



continue evolving as independent Brownian motions after time Ti until the first time 
T2 > Ti when one of them hits the boundary of D. Again at the time T2 the particle which 
approaches the boundary jumps to the current location of a particle chosen uniformly at 
random from amongst the ones strictly inside D. The subsequent evolution of X proceeds 
in the same way. We will say that X constructed above is driven by Brownian motion. 
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The main results in this paper are concerned with Fleming- Viot particle systems driven 
by other processes. 

The above recipe defines the process only for t < t^o, where 



There is no natural way to define the process for t > Too, and, therefore, it is of 
interest to investigate what conditions ensure that Tqo = c>o. In Theorem 1.1 of [5] the 
authors claim that in every domain D G M.'^ and every > 2, we have Too = oo, so 
the Fleming- Viot process is always well-defined. However, the proof of Theorem 1.1 in 
[5] contains an error which is irreparable in the following sense. That proof is based on 
only two properties of Brownian motion — the strong Markov property and the fact the 
the hitting time distribution of a compact set has no atoms (assuming that the starting 
point lies outside the set). Hence, if some version of that argument were true, it would 
apply to almost all non-trivial examples of Markov processes with continuous time, and in 
particular to all diffusions. However, in p], the authors provided an example of a diffusion 
X on D = (0, oo) (a Bessel process with dimension z/ = —4), such that Too < oo for the 
Fleming- Viot process driven by this diffusion with N = 2. 

It is not known whether Theorem 1.1 of [5j is correct in full generality. It was proved 
in [31, [To] that the theorem holds in domains which do not have thin channels. 

1.1. Main results. We will prove two theorems. The first theorem is concerned with 
Bessel processes but it is motivated by the original model based on Brownian motion in 
an open bounded subset of M*^. Recall that for any real u, a //-dimensional Bessel process 
on (0, oo) killed at may be defined as a solution to the stochastic differential equation 



Tr 



OO 



lim Tfc. 



dWt + 



u-1 



dt. 



(1.1) 



t — 



2Xi 



where W is the standard Brownian motion. To make a link between Brownian motion in 
a domain and Bessel processes, we recall that there exists a regularized version p of the 
distance function ([HI Theorem 2, p. 171]). More precisely, there exist < Ci, C2, C3, C4 < 
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oo and a function p : D — t- (0, oo) with the following properties, 

Ci dist(x, < p{x) < C2 dist(a;, 



sup |Vp(x)| < C3, 



sup 



< C4 for 1 < i,m < d. 



The above estimates and the Ito formula show that if i? = (i?^, . . . , B'^) is a (i-dimensional 
Brownian motion and Zt = p{Bt) then 



k=l 



where the functions afc(-) and 6( ■ ) are bounded. This shows that the dynamics of Z 
resembles that of a Bessel process. Note that if Too < C)0 for the Fleming- Viot process 
driven by Brownian motion in a domain D then the distances of all particles to dD go to 
as 1 1 Too; by Lemma 5.2 of [3]. Hence, it is of some interest to see whether a Fleming- Viot 
process based on a Bessel process can become extinct in a finite time. We have a complete 
answer only for N = 2, i.e., a two-particle process. 

Theorem 1. Let X 6e a Fleming- Viot process with N particles on (0, 00) driven by Bessel 
process of dimension z/ G M. 

(i) If N = 2 then t^q < 00, a.s., if and only if u < 0. 

(a) If Nu > 2 then t^o = 00, a.s. 



Our second main result is also motivated by some results presented in Several 
theorems in [5l |3] are concerned with limits when N 00. To formulate rigorously any 
of these theorems it would suffice that t^q = 00, a.s., for all sufficiently large A^. In other 
words, it is not necessary to know whether Too = 00 for small values of A^. One may 
wonder whether it is necessarily the case that Too = 00 for any Fleming- Viot-type process 
and sufficiently large A^. Our next result shows that once the drift of the diffusion is 
slightly stronger than the drift of any Bessel process then Too < 00 for the Fleming- Viot 
process driven by this diffusion and every N. 

Consider the following SDE for a diffusion on (0, 2], 

f' 1 

Xt = xo + Wt- / —^ds-Lt, t<To, (1.2) 
Jo pXf 
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where xq & [0,2], 13 > 2, W is Brownian motion, Tq is the first hitting time of by X, 
and Lt is the local time of X at 2, i.e., Lt is a CAF of X such that 



Jo 

It is well known that (11. 2p has a unique pathwise solution {X,L) (see, e.g., |2j. Theorem 
1.12.1). We will analyze a Fleming- Viot process on (0,2] driven by the diffusion defined 
in (11.21) . The role of the boundary is played by the point 0, and only this point. In other 
words, the particles jump only when they approach 0. Let P'' denote the distribution of 
the Fleming- Viot particle system starting from Xq = x. 

Theorem 2. Fix any (5 > 2. For every N > 2, the N -particle Fleming-Viot process on 
(0,2] driven by diffusion defined in (11.21) has the property that Too < oo, a.s. Moreover, 



where c\ and C2 depend only on N and (3, and satisfy < ci, C2 < oo. 

Remark 1. (i) If we take /3 = 2 in (II. 2p then the diffusion is a Bessel process (locally near 
0). Hence, we may say that Theorem [2] is concerned with a diffusion with a drift "slightly 
stronger" than the drift of any Bessel process. 

(ii) The theorem still holds if the constant 1//3 in the drift term in (II. 2p is replaced by 
any other positive constant. We chose 1//3 to simplify some formulas in the proof. 

(iii) The diffusion (ll.2p is reflected at 2 so that we can prove the exponential bound 
in (II. 3p . For some Markov processes, the hitting time of a point can be flnite almost 
surely but it may have an inflnite expectation; the hitting time of by one- dimensional 
Brownian motion starting at 1 is a classical example of such situation. The reflection 
is used in (ll.2p to get rid of the effects of excursions of the diffusion far away from the 
boundary at 0. A different example could be constructed based on a diffusion on (0, oo) 
with no reflection but with very strong negative drift far away from 0. 

We end this section with two open problems. 




P"(roo >t)< cie 



t > 0, X G (0,2]^, 



(1.3) 



Problem 1. Find necessary and sufficient conditions, in terms of N and v , for non- 
extinction in a finite time of an X -particle Fleming-Viot process driven by u- dimensional 
Bessel process. 
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Problem 2. Does there exist a Fleming- Viot-type process, not necessarily driven by Brow- 
nian motion, such that Tqo = oo, a.s., for the N -particle system, hut Too < oo with positive 
probability for the {N + 1) -particle system, for some N > 2? 

The rest of the paper contains the proofs of the two main theorems. 

2. Proof of Theorem [T] 

2.1. Bessel processes. We start with a review of some facts about Bessel processes and 
Gamma distributions. Let Zf, t > 0, he a square of Bessel process of dimension i/ G M 
starting at x > 0, (Z ~ BESQ'^(x), for short), i.e., Z is the unique strong solution to 
stochastic differential equation 

dZt = udt + 2^/\Z^\dWt, Zo = x, 

where is a one- dimensional Brownian motion (see fTIi Chapter 11] for the case z/ > 
and [S] for the general case). 

Squares of Bessel processes have the following scaling property: if Zf ~ BESQ'^(x) and 
for some c > and alH > we have Z'^ = cZ^-if, then Z' ~ BESQ''(ca;). 

If Z ~ BESQ'^(x) with x > 0, and Tq denotes the first hitting time of 0, then Tq = oo, 
a.s., if > 2, and To < oo, a.s., if u < 2. Moreover, in the latter case we have 

n i |j. (2,1) 

where G is F (l — |) -distributed random variable |i8| eqn. (15)]. Here and in what follows 
we say that a random variable is r(a)-distributed if it has the density 

faix) = 7-r-Ta;"~"^e~^, x > 0, a > 0, 
T{a) 

where 

POO 

T{a) = / x"-^ e-^ dx 
Jo 

denotes the standard gamma function. Note that we consider only a one-parameter family 
of gamma densities, unlike the traditional two-parameter family. 

In [T2j, Bessel process X of dimension u > starting at x > (X ~ Bes'^(x)), is 
defined as the square root of BESQ'^(x^) process Z. If z/ > 0, then by so called comparison 
theorems, the paths of Zt are defined for all t > 0, so Xf is well defined for all t G [0, oo). 
We define Bessel process X of dimension z/ < starting at x > as the square root of 
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a BESQ'^(x^) process Z, i.e., Xf = \fZ't for t < Tq. For any real z/, these definitions are 
equivalent to the definition given in (11. ip by the Ito formula. 

Processes Bes'^(x) with z/ G M scale as follows. If X ~ Bes''(x) is a Bessel process on 
[0,To), then for all c > 0, cXc-2( is a Bes^(cx) process on [0,c^To]. This follows easily 
from the scaling property of BESQ''(x) processes. 

2.2. Proof of Theorem [T] (i). We start with an alternative construction of the Fleming- 
Viot process X. Let X = {Xt, t e [0,To)) be a Bes''(l) process. Let Y = (Y^^,Y^^), where 
and Y^^ are independent copies of Xt and let = i = 1,2,..., be a 

sequence of independent copies of Y. For i = 1,2, . . . we set 

ai = inf {t > : F/'^ A F/'^ = O} , 

and 

It is easily seen that ai may be represented as cJi = min(To,TQ), where Tq is an indepen- 
dent copy of Tq, and that (cTj, z = 1, 2, ...) is a sequence of independent and identically 
distributed random variables. 

We construct a two-particle Fleming- Viot type process = {X^,Xf) as follows. First 
let Ti = 0"! and set X^ = Y^^ for t G [0, Ti). At Ti one of the particles hits the boundary of 
D = (0, oo), and it jumps to = ai. To continue the process we use the scaling property 
of Y^: let T2 = Ti + ^^0-2 and set X^ = ^iY^_2, , for t G [ti,T2). At T2, one of the 
particles hits the boundary and jumps, this time to ^2 = C(2^i- We continue the process 
in the same way by setting 

2=1 

n 

and 

= ^ny%,^., for t G r„+i) . 
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It is easy to see that the construction of X given above is equivalent to that given in 
the Introduction, except that the driving process is a z/-dimensional Bessel process. The 
process is well defined up until Too, and we will show now that Too < oo almost surely 
if and only if z/ < 0. 

Note that Xq = (1, 1) for the process constructed above. However, it is easy to see that 
for any two starting points Xq = (xq,Xq) and Xq = {zq,Zq) with xl,xl, z^, Zq > 0, the 
distributions of X^^ are mutually absolutely continuous. This implies that the argument 
given below proves the theorem for any initial value of X. 

The case > 2 is very simple: then ai = oo, a.s., so t^c = oo, a.s. So for the rest of 
this section we assume that u < 2. 

To check whether Too < oo or Too = oo, we will apply the following theorem. Let 
log^ X = max(logx, 0). 

Theorem 3. (\Q\; see also P] or [9\) Let n>l} he a sequence of independent 

and identically distributed random vectors such that G M and 

E (log+ |Ai|) < oo, E (log+ < oo. 

Then the infinite random series 

oo n— 1 
n=l j=l 

converges a.s. to a finite limit if and only if 

E log I All < 0. 

We will apply Theorem [3] with An = ot^ and -B„ = cr„ . Thus, in order to prove Theorem 
d] (i), it suffices to show that 

(i) E log (Ti < oo for z/ < 2; 

(n) Elog(af) < for z/ < and E log(a?) > for z/ > 0. 
Proof of (i). Note that, in view of fl2.ip . 

Elogdi < ElogTo = Elog-^ = -log2 -ElogG, 
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where G ~ F (l — l) . But for G ~ r(a) with a = 1 — | > 0, we have 

/>oo 

ElogG= / \ogx fa{x)dx 
Jo 



oo 



r(a) Jo 

1 d 



T(a) = ^^(a) < oo, 



r(a) da 

where il) is well known digamma function ([1] Section 6.3]) defined as 

V'(a^) = -J- logr(x) ^ 



dx r(a:) 

Proof of (a). By Theorem 11 of [TT] we get that the density of a\ is given by 

(y + 2r3 - r(3-z. + 2n) / y y 
^ r(i-i) A.,,r(2-| + n) l(y + 2)2j 

r(i-|)n(i/ + 2)^ 



where 



with 



oo 

.Z" 

n=0 



r (2n + 3 - I/) 
= -^^ -_ ^, n = 0,l,2,.... 



^" n!r + 2 - I) 
By the duplication formula for the gamma function ([H eqn. 6.1.18]), i.e.. 



we have 

_ 22^^+2-^ r (n + ^) 



2 



r 



,'7r \ n 
where for a; > —1 and G N, 

^x\ r(x + 1) 



^fcy A;!r(x-A; + 1) 
is a generalized binomial coefficient. Therefore, 

n=0 



TT V 2 / Vl-42 
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as for a G M 



n=0 

Now 



(t, + 2)'V ^ \ 2 ) (j,2 + 4)^ 



and therefore for y > 

r(^) 



v^r(i-i) (^2 + 4)^- 

So, to prove (ii) we need to study the sign of the integral 

/(z/) = / K{y) logydy. 
Jo 

Recall the Student's t-distribution with a > degrees of freedom (jTJ section 26.7]). The 
density for this distribution is given by 



^ \ 2 



fix; a) = ^ 2 / 1 H 1 , -oo < a; < cx). 

' ^ v^r(f) \ a ' 



Changing the variable y = ,T^_ in /(z/) we get 



2-u 

2x 



poo 

I{u)= / /(x;2-z/)log 
Jo 



: C/X 



1^, 21X1 
-E log ■ 



2 ° v/2^' 

where X is a random variable with t-distribution with (2 — i^)-degrees of freedom. 
It is well known ([li section 26.7]) that 



Jv — — 



V ' 

where Z has standard normal distribution and V has chi-squared distribution with (2 — v) 

degrees of freedom, and Z and V are independent. Therefore 

1 2 \Z\ 
= -Elog' 



2 " 

1 / V 
log 2 + - E log — - E log — 



Note that ^ has T (i) distribution and ^ has F (^-g^) distribution. Therefore, by fl2.2p . 

'2 - z/^ 



/(z/) = -log2 + -U'U 
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The function ifj is strictly increasing with ip = —2 log 2 — 7 and ^^(l) = —7 where 7 is 
the Euler constant (PQ eqns. 6.3.2, 6.3.3]). Using these facts we see that 



4 V V 2 

and therefore liy) < iff > 1 iff z/ < 0. This completes the proof of (ii) and of 
Theorem [1] (i). □ 

2.3. Proof of Theorem [1] (ii). Suppose that X = (X^,...,X^) is a Fleming- Viot 
process driven by i/-dimensional Bessel process, Xq = (x^, . . . , x^), > for all 1 < 
J < N, and Nu > 2. Let Zt = (X/)^ + ■ ■ ■ + {X^f and zq = {x^Y + ■■■ + (x^)^ > 0. 
According to [131 Thm. 2.1], the process {Zt,t G [0,ri)} is an (A^z/)-dimensional square 
of Bessel process, i.e., it has distribution BESQ^^(2;o). More generally, {Zt,t G [rfc,rfc+i)} 
has distribution J B'ESCl^'^ {z)¥(Zr^ E dz) for A; > 0, where, by convention, tq = 0. Let 
Yt = Z^ , Bq = and define B inductively on intervals (r^, r^+i] by 

B, = B,,+ dY, - / ^^ds. 

Then, by the Ito formula, i? is a Brownian motion and 

dZt = 2^tdBt + Nvdt, 

for t G (rfc, Tfc+i), /c > 0. Let be defined by Zq = zq and 

Zt= I 2\ Zs dBs + Nut, t > 0. 





By definition, Z is an (A^z/)-dimensional squared Bessel process on [0,oo). We assumed 
that Nu > 2 and 2:0 > so we have Zt > for all t > 0, a.s. Since Z is continuous, 
for every integer j > there exists a random variable aj such that Zf > aj > for all 
t G [0,j], a.s. Note that Zt = Zt for t G [0,ri) and Z^-^ < Z^-^ because Z has a positive 
jump at time ti. Strong existence and uniqueness for SDK's with smooth coefficients 
implies that Zt < Zt for all t G [ti,T2), because if the trajectories of Z and Z ever meet 
then they have to be identical after that time up to T2. Once again, Z^.^ < Z^.^ because 
Z has a positive jump at time T2. By induction, Zt < Zt for all t G [Tk,Tk+i), k >0, a.s. 
Hence, Zt > aj > for all t G [0, j] and j > 0, a.s. This implies that Too = 00, a.s., by an 
argument similar to that in Lemma 5.2 of [5]. 
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3. Proof of Theorem [2] 

3.1. Preliminaries. We will give new meanings to some symbols used in the previous 
section. Constants denoted by c with subscripts will be tacitly assumed to be strictly 
positive and finite; in addition, they may be assumed to satisfy some other conditions. 

(i) Let W be one- dimensional Brownian motion and let 6 be a Lipschitz function defined 
on an interval in M, i.e., \b{xi) — b{x2)\ < L\xi — X2I for some L < 00 and all xi and X2 
in the domain of b. Consider a diffusion Xt,t G [s,u], satisfying the following stochastic 
differential equation, 

dXt = dWt + b {Xt) dt, Xs = a. (3.1) 
Let i/t be the solution to the ordinary differential equation 

^yt = b{yt), Vs = a. 

We will later write y' = b{y) instead of ^yt = b{yt). 

The following inequality appears in Ch. 3, Sect. 1 of the book by Freidlin and Wentzell 
[7]. For every 6 > 0, 

P ( sup \Xt ~yt\>s) < P ( sup \Wt\ > Se'^^^'-A , 

\s<t<u J \s<t<u J 

where L is a Lipschitz constant of b. It follows that 

P ( sup \Xt-yt\>5] <¥{ sup Wt > 5e"^(""'M + P ( inf Wt < -5e~^("""M 

\s<t<u J \s<t<u J \s<t<u J 

= 2P ( sup Wt > 5e-^("-^M 

\s<t<u J 

= 4P (W„ -Ws> (5e-^("-')) 

where Cq is an absolute constant. 

(ii) Recall that /3 > 2 and consider the function 

&(x) = --^, x>0. (3.3) 

We need the assumption that /3 > 2 for the main part of the argument but many calcu- 
lations given below hold for a larger family of /3's. It is easy to check that 

ys,a{t) := (a^ + s - t) s<t<s + a^, (3.4) 
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is the solution to the ordinary differential equation 

V = b{y) (3.5) 

with the initial condition ys,a{s) = a, where s G M, a > 0. Note that the function ys,a{t) 
approaches vertically at t = s + a^. 

(iii) Fix any 7 G (0, 1) and let L be the Lipschitz constant of b on the interval 
[a(7/2)i/^/2,2a]. Then 

L = 6'(a(7/2)V/^/2), = 

and, therefore, 

(3.6) 



/3(a(7/2)V/3/2)/3 

Let X be the solution to (13. ip with h defined in (13. 3p . Assume that 5 > is so small that 

0(7/2)^/^^/2 < yo,a((l - 7/2)a^) - 5 < 2/o,a(0) + 5<2a. (3.7) 

It follows that if supo<t<(i_^/2)a/3 \Xt - ya,ait)\ < 5 then Xt G [0(7/2)^/^/2, 2a] for < 
t < (1 — 7/2)a^. Hence, we can apply (13. 2p with L given by (13. 6p to obtain the following 
estimate 

P( sup |Xt-yo,a(t)| >5|Xo = a| <Coexpf-cAy (3.8) 



0<t<(l-7/2)a/3 



where Xt satisfies (13. ip and Ci depends on {3 and 7, but it does not depend on 5 and a. 
(iv) Suppose that a > and u = (1 — 7)0''. Then yo^a{u) = 07^/''. For e G (0, 1), let 

6 = 6{e) = a[l-{l-e^y/^]. (3.9) 

We fix e G (0, 1) so small that for all a > the inequality (13. 7p is satisfied with S in place 
of 6 and, moreover, 

(l-7/2)(a-5)'' > (l-7)a^ (3.10) 

and 

7^/^ + 1 - (1 - £7)1/^ < 1. (3.11) 

If (5 G [-5, 6] then 

yo,a+s{t) = ((a + 5)'' - t) < t < (a + S)f'. 
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It is straightforward to check that 



We will estimate the difference yo^a+siu) — yo^a{u) as a function of 6. Let 

f{S) = yo,a+5{u) = ((a + - (1 - 7)a^) • 

Then 

yo,a+siu) - yo,a{u) = f{5) " /(O) = 5f'{5), 
where S is between and S. But 

f'{S) = {a + 6)'-' ((a + 6)' - (1 - 7)a^)'^'~' 

and 

r (5) = -a^(/3 - 1)(1 - 7) (a + 6^'' ((a + 6)^ - (1 - 7)^^) 

It follows from (^M) and the fact that e E (0, 1) that (a + 5)^ - (1 - 7)a^ > 0. Thus 
/" < and /' is strictly decreasing. Therefore, for 6 G [—6, 6] we have 

/'W</M)=^r^ =:M>1. (3.12) 



1//3-2 



.7(1 

Hence \yo^a+5{u) — yo,a 

(m)| < M\6\ for 5 G [-5,5]. Suppose that 6' G (0,5] and 6 - ^ 



A/+1- 

If ao G [a — 5, a + 5] , then 

yoA^) - < l/0,ao(M) < Z/0,a(M) + 5'. (3.13) 

The function b{x) is strictly increasing on (0, 00). This easily implies that if < ai < 02 
and < s < t < then 

2/0,a2(s) -?/0,ai(s) < Z/O.aaW -2/0,ai(^)- (3.14) 

Hence, inequality fl3.13p holds in fact for all t G [0,u] in place of u and, moreover, 
yo,a-s{'t) < yo,a(t) — 5'. Another consequence of (13.131) and (I3.14p is that if s G [0,u] 
and ao G [l/o,a(s) - 5,?/o,a(s) + 6] then for t G [s,m], 

l/0,a-5W < 2/0,a(t) -6'< ys,a,{t) < l/o,a(t) + 6'. (3.15) 
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The following generalization of fl3.10p also follows from fl3.14p . If s G [0,u] and oq > 
l/o,a(s) - 6 then 

s + (l-7/2)a^ > (l-7)a^. (3.16) 

3.2. Proof of Theorem [21. Suppose that = (X/, . . . , X/^) is a Fleming- Viot process 
on (0,2] driven by the diffusion defined in (11. 2p . with an arbitrary 2 < N < oo. Recall 
that the role of the boundary is played by the point 0, and only this point. In other 
words, the particles jump only when they approach 0. 

Step 1. Let X = (x^, . . . , x^) , [N] = {1, . . . , N}, and let ji be the smallest integer in 
[N] with 

x-'^ = max X-' . 
i<i<JV 

Consider process X starting from Xq = x and let Ii = {ji}, Ji = [N] \ Ii and 5*0 = 0. 

Let u = (1 -7)(x-'i)^ and 

^1 = M A inf {t > : 3j^j^ X/ = X/^ } . 

Note that two processes X' and can meet either when their paths intersect at a time 
when both processes are continuous or when one of the processes jumps onto the other. 
Let j2 be the smallest index in Ji such that the equality in the definition of 5*1 holds with 
j = ja- Let h = {ji, ja} and Ja = [N] \ h- 

Next we proceed by induction. Assume that, for some n < N, the sets 
Ji, . . . , Jn, and stopping times < ^2 < . . . < 5'„_i are defined. Then we let 

S^ = uA inf {t > : 3,g,„ 3,gj„ X^ = X/} , 

In+i = /n U {jn+i} and Jn+1 = [N] \ In+1, where jn+i is the smallest index in J„ such that 
the equality in the definition of Sn holds with j = jn+i- 

The set has n elements which are indices of particles which are "descendants" of the 
particle X-'^ that was the highest at time 0. By convention, we let /„ = In and Sn = u 
for n> N. 

Step 2. Write a = x^^ and u = (1 — 7)0^. Then x G (0, a]^. Recall 6 and M defined in 
([33]) and ([312]) . and for 1 < n < X define 
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Note that 5„ = (M + Consider events 



F„ = y i sup \xi - ?/o,a(t)| > L \ 

■ J \^Sn-l<t<Sn ) 



Note that for every t, maxj-g/^ X/ > maxjgj,^ Xf. Hence, 



U I sup X{ - yoAt) >Ac [j Fn, 

l<j<N '^°^*<" ^ l<n<N 



and, therefore. 



P'^ ( U I sup Xi - yoAt) > 4 

\l<n<Af / 

= P"I U FnnF^n---nF^_, 

\l<n<N , 

< J2 (i^n n n ■ ■ ■ n i^^i) 

l<n<Af 

< Yl p^(i^n|F^,...,F„%) 



(3.18) 



l<n<Af 



^.^„<' M AlzT \Sn-\<t<Sn 



Fl, . . . , -F^„l , 



l<n<Af jel„ 

where we adopted the convention F^{Fi \ Fq) = F^{Fi). 

Suppose that F^_^ holds and j G /„. Then \Xg^__^—yo^a{Sn-i)\ < ^n-i- Let yl,t > Sn-i, 
be a solution to y' = b{y) with ys^_^ = -^Sn-i- (13.151) . 



P'' sup XI - yoAt) >Sn\ F:_i < PM sup Xf - yl > 5„ | F^^ . 

\S^-l<t<Sn J \S„-l<t<Sn / 

It follows from (13.161) that we can apply (13.81) (with an appropriate shift of the time scale) 
to X^ , assuming that holds, on the interval [^„_i, Sn] C [0,u] = [0, (1 - 7)0^]. We 
obtain 

P'' ( sup \Xi-yl\>6n\ F„=_i ) < Co exp I -c^ 



5„_i<i<5„ / y (2/0,a(5'„-l) + 5„-l)^ 

< Co exp -C2 — — = Co exp -C3 



a + 5y J "V (a + S) 
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We combine this estimate with f l3.18p to see that 



f 1) j sup 



Xi - yoAt) > S 



(3.19) 



\Sn-l<t<S n 



l<n<N jeln 



< CqN'^ exp -C3 



52 



Step 3. We will prove that there exist v < 00 and r G (0,2) such that if x G (0,r]^ 
then 



P"(r^ >v)< 



1 



(3.20) 



Consider an r G (0,2) and for x = (x^, . . . ,x ), let = max^x-^, Uo = 0, and for 
A; = 0, 1,2, let 



Uk+1 = f/fc + (1 - 7)^: 



Ak+i = max Xjr 

Let Y^'^ denote the solution to ODE fl3.5p with the initial condition Y^^ = A^. Recall e 
from flXTUl) and let = Ak[l-{1- £7)^/^] . For A; = 0, 1, 2, . . . define events 

max sup xi >Y^ + A 

l<i<A'C7fc<t<{/fc+i 

Note that Y^}^^^ = -f^^^Ak. Suppose that 0^=0 ^fc holds. Then 

Ak+i < 1^1^ Ak + /I41 - (1 - £7)1/^] = C4Afc, 

for all fc, where C4 = 7^/^ + 1 - (1 - £7)^/^ < 1, by fIXTT]) . Hence, < c^Aq and, 
therefore, X^fc^f ^ -'■^ f = (1 — 7) ^^li '"''£4'' < 00 then limfc_^oo f^/t < 

(1 — 7) X^fcLi ^Q^l — ^ ^^"^ limsup^^^ maxi<j<jvX/ = 0. This implies easily that Too < v. 
Thus, to prove (13.201) . it will suffice to show that there exists r G (0,2) such that if 
X G (0,r]^, then 

(3.21) 



\k=0 / 



But 



00 

p-( IJrJ <p-(ro) + 5^P''(r, |rs,...,p,_J. 



fe=i 



(3.22) 
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By fl3.19p and the strong Markov property applied at f/^, 

r(r„|rs.....n_0<-o'v=exp(-c,^^^-^) 

<coN'exp[-c,Al-^c1^'-'^ 
where C4 < 1. So by (13.221) . if max^ < r, then 



5 



\k=0 J k=0 

which is convergent. Since 2 — /3 < 0, we can choose r > so small that the above sum 
is less than 1/2, proving fl3.2ip . 

Step 4- Let r G (0,2) and v be as in Step 3. Partition the set (0,2]^ into two sets 
A = (0,r]^ and A'^. First we will show that the time when process X enters the set A 
has a distribution with an exponentially decreasing tail. 

So assume that Xq G A'^ and let 

h = {je[N]:X',eir,2]}, h = [N]\h. 
Let t( be the the first hitting time of by the process and let 

if /2 

maxjg/2 {'''1} ' otherwise 
Consider 

pi(x) = {V,e/, Vo<*<i/2 Xi G 2] ; < 1/2; Wi^iXl<t<i/2 e 2] } . 

We will argue that for x G we have 

Pi(x)>pi>0. (3.23) 

Indeed, with probability at least gi > any particle from Ji stays in the interval [r/2,2] 
up to time t = 1/2. With probability at least ^2 > any particle from I2 hits before 
time t = 1/2; with probability at least 1/A^ it jumps onto a particle in Ji; and then with 
probability at least ^3 > it stays in the interval [r/4,2] up to time t = 1/2. Therefore 
(I3.23P holds with pi = {qiq2(l3/N)^ . Obviously qi, q2, qs and pi depend on r. 
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Next, if we define 

and B = [|,2] , then an argument similar to that proving f l3.23p shows that for ^ E B 
we have P2(x) > P2 > 0, where p2 depends on r. Therefore, by the Markov property at 
time t = 1/2, for x G we have 

P^(Xi eA)>p:= P1P2 > 0. (3.24) 

Now let 

T = inf {t > : G A} . 

By ([321, for all x G (0,2]^, 

P''(T < 1) >p > 0. 
Applying the Markov property at t = 1, 2, . . . we obtain 

P^(r > A;) < {l-pf. 

Choose k so large that (1 —p)^ < ^. Recall that r and v are as in Step 3. Let 9 denote 
the usual Markovian shift operator. Then for any x G (0, 2]^, 

P'' {Too>k + v)< P^(T > A;) + P^(roo o Or > v) < {1 - p)'' + - := q < 1. 
Therefore, applying the Markov property at times k + v, 2{k + v), 3{k + f ),..., we obtain, 

P"(roo >r2(A; + t;)) <g", n = l,2,..., 
which proves (11.31) . This implies that Too < cxd, a.s. 
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